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ABSTRACT 



A horizontal cylinder, located near the floor of a two 
dimensional wave channel was subjected to a train of non- 
linear gravity waves. The horizontal and vertical forces 
were measured and presented in dimensionless form. Experi- 
mental values of horizontal and vertical force coefficients 
are presented as functions of dimensionless wave height and 
dimensionless wave period. The dimensionless force coeffi- 
cients predicted by a modified Morrison equation and a 
Froude-Krilov force are compared to experimental data. 

Fluid particle velocity and acceleration values were calcu- 
lated from Stokes fifth-order gravity wave theory. Experi- 
mental dimensionless wave periods from 10 to 200 were 
investigated. 

The horizontal force coefficients were found to vary 
linearly with wave height for dimensionless period values 
from 6 0 to 120. The vertical force coefficients were 
found to be inertia dominated at low dimensionless wave 
heights and dominated by a lift force at higher wave heights. 
The theory predicted experimental values of force coefficients 
with good accuracy, especially at greater water depths. 
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I. INTRODUCTION 



Exploration for and use of ocean resources has received 
increasing attention in the scientific community in recent 
years. A considerable effort has resulted in the development 
of technology in areas such as off-shore oil exploration, 
recovery of petroleum and deployment of underwater habitats. 
Oil production often involved the use of very large submerged 
oil storage tanks and associated piping systems. In some 
locations, oil is transferred long distances in submerged 
piping systems, often leading through the surf zone to 
facilities ashore. Also, waste disposal systems usually 
involve deployment of large diameter outfall pipes which are 
laid through the surf zone. Proposals for large scale mining 
and food gathering activities indicate the trend toward 
design and construction of large submerged structures. 

This activity has generated considerable interest in the 
interaction of gravity waves with submerged objects and 
particularly in the forces induced. Most previous investi- 
gations [1,2,3] in this area center on the use of a so-called 
Morrison Equation [3], which involves both a drag and iner- 
tial component of force. Unfortunately, the drag and added 
mass coefficients are found to be neither constant nor 
simply shape - or Reynolds number - dependent in the unsteady 
flow that occurs in the interaction of waves with fixed 
bodies. In addition, where the amplitude of wave motion is 
large compared to the body dimensions, flow separation occurs. 
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and the coefficients vary as the flow in a given direction 
develops. Finally, for small depth to characteristic dimen- 
sion ratios, the proximity of the free surface may influence 
the value of the coefficients, even if no flow separation 
occurs . 

In spite of these difficulties, the Morrison Equation 
will, with certain restrictions, provide a complete descrip- 
tion of the horizontal forces acting on a large horizontal 
cylinder in contact with a plane boundary. If the cylinder 
is large in comparison to the amplitude of the fluid motion, 
separation does not occur, and the resulting wave-induced 
oscillatory flow may be considered an unseparated potential 
flow about the cylinder. This has been verified by Kuelegan 
and Carpenter [2] and by Sarpkaya and Garrison [4], For 
such cases, the drag component of force is negligible, while 
the added mass coefficient is a constant. 

Application of the Morrison Equation to determine the 
vertical components of force for a bottom-mounted cylinder 
yields a vertical force component of zero, an obviously 
invalid result. It is recognized that the vertical force 
should be composed of two components; one associated with 
the inertia effects and a second which accounts for the lift 
force caused by the increased velocity and hence decreased 
the pressure over the top of the cylinder. 

The inertial component of vertical force can be approxi- 
mated by what is commonly known as a Froude-Krilov force. 
This is the force caused by the pressure distribution around 
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the surface of the cylinder existing in the wave if the 
cylinder was not present. This is an approximation which 
does not account for the presence of the body and is 
therefore considered to be an underestimate of the force. 

The lift component of vertical force can be approximated 
by use of the unseparated, potential flow model used in the 
analysis of the horizontal forces, i.e., that of a uniform 
flow past a cylinder in contact with a rigid boundary. The 
lift coefficient for this case has been determined by Dalton 
and Helfenstein [5]. 

A study carried out by Johnson [6] for horizontal forces 
acting on a bottom-mounted, horizontal cylinder in long waves, 
assumed that at long wave lengths, the horizontal forces 
depended on wave height and water depth only. This assumption 
is, however, of doubtful validity in the range of Johnson’s 
test . 

Shiller [ 7 ] measured wave forces on a submerged horizontal 
cylinder due to small amplitude waves. Over the range of 
wave heights considered, the magnitude of the horizontal 
force was found to vary linearly with wave height. The 
magnitude of the vertical force was found to increase in 
proportion to the wave height squared. 

Perkinson [8] extended Shiller* s investigation to include 
larger wave heights and periods. His study indicates that 
horizontal forces increased linearly with wave heights and 
becomes independent of wave period for large wave periods. 

The vertical forces were found to contain two regimes: the 
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lower wave lengths where lift force predominates; and the 
higher wave lengths where inertial force is dominant. 

In order to evaluate any forces using the above analyses, 
it is necessary to first evaluate fluid velocities, accelera- 
tions and pressure by use of some wave theory. Fairly simple 
expressions may be derived for the forces by use of Airy wave 
theory. However, for finite amplitude waves in finite depth 
water, a non-linear wave theory is considered to be necessary. 
Stokes fifth-order wave theory has been chosen for this 
study. 

Stokes originally developed a second-order theory [9] for 
the case of a non-linear wave in water of finite depth. This 
method has been extended by Borgman and Chappelear [10] to 
third order. Skjelbreia and Hendrickson [11] have extended 
the solution to third and fifth order. Bretschneider [12] 
has presented a method for extension to any order. 

The purpose of this study is to compare the analytical 
results of Stoke* s fifth-order wave theory with experimental 
values at longer wave periods. The transition section of the 
wave channel used by Perkinson was extended approximately 
forty feet to provide a slope of 1:20. This provided the 
possibility of extending Perkinson’ s work to include longer 
wave periods. This study, then, is an extension of the 
work of Shiller [7] and Perkinson [8], 
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II. THEORETICAL ANALYSIS 



A. PROBLEM DEFINITION 

The problem under consideration is illustrated in Figure 
1. A train of gravity waves is considered to progress in 
the positive x-direction in water of depth h, with rate of 
propagation C. The fluid particle velocity is expressed in 
terms of the horizontal and vertical components, u and v. It 
is of primary interest to determine the horizontal and 
vertical components of wave force acting on the horizontal 
cylinder which is in contact with the rigid bottom. 

B. DIMENSIONAL ANALYSIS 

An exact analytical solution to this problem is quite 
formidable and, therefore, an approximate solution is con- 
sidered in this work. However, first it is instructive to 
carry out a dimensional analysis of all pertinent parameters 
involved. 

In general, the maximum of the wave force per unit 
length acting on a submerged cylinder due to wave motion is 
known a priori to depend on the following parameters: 

P max A = ^(h.H.T.L.a.p.g.y) (D 

where 

F = wave force 
h = water depth 
H = wave height 
T = wave period 
L = wave length 
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T = Wave Period 




Fig. I. DEFINITION OF GEOMETRY 



14 



a = cylinder radius 
p = fluid density 
g = gravitational constant 
y = fluid viscosity 
i = cylinder length 



However, a relationship exists between the parameters 
associated with the incident wave (i.e. h,H,T and L). 
Consequently, only three of these parameters (e.g. h,H and 
T) are needed to completely specify the incident wave. A 
dimensional analysis of the parameters in equation (1) yields 
the following groups: 

F ma x // PS a2jl= i ' 2 ^ sT2//h> h//a > H//2a > y/^"^gha 2 ) (2) 



The last term on the right side of equation (2) repre- 
sents the ratio of Proude number to Reynolds number, indi- 
cating the ratio of viscous to inertial forces. It is 
believed that, if this number is small, it may be neglected 
from further consideration, and equation (2) may be rewritten 



max 



= f ^ (gT /h , H, h) 



where 



max 



H 

h 



P /pga £ 
max 

H/2a 

h/a 



(3) 



The approximate analytical approach of this study seeks: 
1. To develop expressions for the fluid particle motion 
in the gravity wave as functions of the incident wave 
parameters, and 
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2). To develop expressions for the forces on the sub- 
merged cylinder as functions of the fluid particle 
motion. 

The assumptions upon which this analysis is based are: 

1) . Water will be considered invlscid and incompressible, 

2) . The cylinder radius will be considered small in 

comparison to the dimensions (h and L) of the 
incident wave. 

C . WAVE THEORY 

The motion of fluid particles in a gravity wave based on 
assumptions 1) and 2) above is specified by the following 
boundary value problem. 

1). Governing Equation 

Assumption 1) implies that Laplace's Equation 
governs the fluid motion. That is: 

V 2 <1> = 0 (4a) 



where <j> is the velocity potential and the two velocity 
components are given in terms of <J) as 



u 



8x * 




(4b) 



2). Boundary Conditions 

The boundary condition on the bottom is that of a 
non-porous wall, i.e., 



8 (j) 

3y 



0 at y = -h 



(5) 
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Two boundary conditions are encountered at the water 
surface, one dynamic and one kinematic. 

The dynamic boundary condition is obtained from 
Bernoulli’s Equation by setting the pressure equal to zero. 
That is: 

u 2 + v 2 + 2-||- = -2g(n - y) at y = n (6) 

where II is the total energy head, a constant. 

Assuming the wave travels without change in form, it is 
then possible to choose a reference system moving in the 
positive x direction with wave celerity C. This makes the 
fluid motion steady with respect to the moving reference 
system. The Bernoulli Equation in such a reference system 
is : 

(u - C) 2 + v 2 = 2g(n - y) at y = n (7) 

The kinetic boundary condition imposes the condition that 
no fluid be transported across the free surface. Or, in 
other words, the free surface must be a streamline. This may 
be expressed mathematically in the form: 

at y = n (8) 

Equations (4a), (5), (7), and (8) completely specify the 
boundary value problem. 

There are, historically, several approximate solutions 
to this problem. The solution to be considered further in 
this analysis is that commonly known as Stokes fifth-order 
wave theory. 
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D. STOKES FIFTH ORDER SOLUTION 



Stokes proposed [9] a power series type of solution to 
the above problem and proved the validity of the second order 
solution based on this power series. Skjelbreia and 
Hendrickson [9] have extended this theory to fifth-order, 
yielding a solution to the boundary value problem of the 
form: 

4> = ^-{A.. cosh($S)sin0 + A o cosh(2gS)sin(20) 

+ A2COsh(33S)sin( 30 ) + A^cosh(43S)sin(^0) 

+ A 5 cosh(5eS)sin(59)} (9a) 

where A^ through A,- are as given in Appendix A. 

This equation may be written in dimensionless form as: 

<ji = C{A, cosh(£S)sin0 + A_cosh(23S)sin(20)+. . . } 

1 2 (9b) 

where: 4> = 3(f>/"\/gh 

C = c/yih 

Differentiation of equation (9a) with respect to time 
yields 

H = C 2 {A cosh(BS)cos0 + 2A 2 cosh(26S)cos(20) 

+ 3A^cosh(3BS)cos ( 36 ) + 4A2|Cosh( i igS)cos(40) 

+ 5A^cosh(53S)cos (56 ) } (10a) 
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or , in dimensionless form 



!£ = C 2 {A 1 cosh(6S)cos(9) + 2A 2 cosh(2gS)cos(20) 

+ . . . } (10b) 

where if = ft / '/sh 

The series form of the wave profile is assumed to be: 

y = ^-{Xcos(0) + B 2 cos(20) + B^costSO) 

+ B^cosClJO) + B^cos(59) } (11a) 

where B 2 through B^ are as given in Appendix A and X is an 
arbitrary constant. 

Equation (11a) can be expressed in dimensionless form 
as 

y = —{X cos(0) + B~cos(20) + ...} (lib) 

d 

where y = y/d 

d ** 2-rrh/L 

The series form for n is assumed to be: 

n = i{x 2 c 3 + X 4 C 4 } (12) 

where and are as given in Appendix A. 

The wave celerity expressed in series form is 

C = | C 2 {1 + X 2 C 1 + X4C 2 > (13a) 
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or, in dimensionless form 



C = r C 2 {1 + X 2 C, + X4C^ } (13b) 

d ° ± d 

where and C 2 are as given in Appendix A. 

Using equations (4b) and (9b) we find the horizontal 
particle velocity to be: 

u = (f^O = C{A n cosh( £S) cos ( 0 ) + 2A o cosh(23S)cos(20) 

d X A. d 

+ SA^coshOBS^osOB) + 4A^cosh(4BS)cos (40 ) 

+ 5A 5 cosh(53S)cos(50)> (14) 

where u = u/y/gh* 

Also, from equations (4b) and (9a): 

v = (||-) = C{A 1 sinh(eS)cos(0) + 2A 2 sinh(2gS)cos(20) 

+ 3A^sinh (33S)cos(30) + 4A^sinh(43S)cos(40) 

+ 5A 5 sinh(53S)cos(50)} (15) 

where v = v/\/gh 

Differentiation of equations (l4) and (15) with respect 
to time yields the particle accelerations 

a = (ttt-) = dC 2 {A,sinh(0S)sin(0) + 2A o sinh(2$S)sin(20) 

X d o _L d 

+ 3A2sinh(3BS)sin(30) + 4AjjSinh(4BS)sin(40) 

+ 5A 5 sinh(53S)sin(50) (16) 

where a x = a x /g 
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and 



= dC 2 {A^sinh(3S)cos(0) + 2A 2 sinh(2£S)cos(20) 

+ ...} (17) 

where a = a /g 

The potential function and it’s derivatives, as well as 
the wave profile are now known in terms of the unknown 
quantities X and d. Specifying equations for these parameters 
are found by developing equations for the peak-to-trough 
wave height and the wave period. That is, the dimensionless 
wave height is defined as : 

fid = J e=0 - y e _„ (18) 

where Hd = H/h 

Substitution of equation (13b) into equation (18 ) 
yields : 

Hd = 2 {X + X 3 B 33 + X 5 (B 35 + B 55 )} (19) 

where B 33 , B^, and B^^ are as given in Appendix A. 

Now, the wave period and the wave celerity are related 
by the relationship: 

C = L/T (20) 



So, 



equation 

(L/T) 2 



(15a) becomes 




p 



( 21 ) 



21 



or ^ 



h/gT 2 = — —75 — C 2 (l + X 2 C, + X H C 0 ) (22) 

4n 2 0 1 2 

2 

The parameter gT /h is referred to herein as the period 
parameter . 

Equations (19) and (22), if solved simultaneously, yield 
- ~ 2 

values for X and d in terms of Hd and gT /h, quantities which 
are incident wave parameters. Thus the potential function, 
particle velocities and accelerations are shown to be functions 
of the dimensionless parameters obtained by dimensional 
analysis and given by equation (3)* (Note that Hd = 2H/h) 

The next step in the analysis is to express the forces 
acting on the cylinder as functions of the potential function, 
particle velocities and accelerations. 

E. FORCES ON THE CYLINDER 

The forces acting on the cylinder are described in terms 
of their horizontal and vertical components. 

1). Horizontal Component 

The horizontal component of force acting on the 
cylinder is expressed in the form of the so-called Morrison 
Equation as . 

C 

F x = "jr( P2a&u 2 ) + (1.0 + C M )Tra 2 *a x (23) 

where C^ = drag coefficient 

C^ = added mass coefficient 
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For cases where fluid particle motion is small in 
comparison with the cylinder diameter, flow separation does 
not occur and the contribution of drag to the total force 
may be disregarded. In this case, equation (23) may be 
rewritten 



The added mass coefficient for a circular cylinder in 
contact with a rigid wall was given in closed form by 
Garrison [13] as = 2.2S. 

2). Vertical Force Component 

As mentioned previously, application of the Morrison 
Equation to the vertical component of force yields a result 
which is clearly erroneous. The vertical force may be 
considered to arise from two sources: an inertial force 
and a lift force. 




(2A) 




a). Lift Force 



The lift force may be expressed as: 



F l = l/2(2a£p)C L u|u| 



(25) 



where 



C 



L 



lift coefficient 



or, in dimensionless form: 



F l = hC L u|u| 



(26) 
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Dalton and Helfstein indicate [7] that C =4.49 



for unseparated flow past a circular cylinder in contact 
with a rigid wall. 



b). Inertial Force 

The pressure distribution around the surface of 



the cylinder existing if the cylinder is not present results 
in a net vertical force (known as a Froude-Krilov force). 



acting on a differential area of the cylinder dA as shown 
in Figure 2. The above equation, written in dimensionless 
form, is: 

_ 2tt _ 

F T = / P(a,y)sinydy (30) 

L 0 

The pressure is obtained from the Bernoulli Equation: 



Consider the resultant differential force dF. 

i 



P = -p{|(u 2 + v 2 ) + |£} 



(31) 



or, in dimensionless form: 



P = -h[|(5 2 + ^ 2 ) + (§£)] 



(32) 



where P = P/pga 



Therefore 




) + (|£)]sinydY 



(33) 



And the total vertical force may be written: 
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Fig. 2 

INERTIAL FORCE GEOMETRY 



I 
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F. DISCUSSION OF COMPUTER PROGRAM 



The equations developed above were programmed for a 

digital computer as given in Appendix D. The basic inputs 

~ p 

to the program are H, h and gT /h. The program yields 
particle velocity, acceleration and computes the horizontal 
and vertical forces on the cylinder at various wave angles. 

The STAKES5 subroutine calculates the coefficients for 
the Stokes fifth order equations. Simultaneous solution 
to equations (19) and (22) are obtained by assuming an initial 
value of X and d based on linear wave theory and iterating 
to the actual value using the Newton-Raphson technique. 

The main program generates tables of values and in-line 
graphs of the wave profile, particle velocity vs. depth, 
vertical and horizontal force vs. phase angle. 
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III. DESCRIPTION OF APPARATUS AND 



EXPERIMENTAL PROCEDURE 



A. WAVE CHANNEL 

Wave forces on the cylinder were measured experimentally 
in a wave channel. The basic wave channel, as shown in 
Figure 3, was fifteen inches wide and consisted of two sec- 
tions: a transition section and a shallow water section. 

The overall channel length was 87 feet. 

Three-quarter inch exterior plywood sheets were used for 
the channel sides. Vertical rigidity was provided by bracing 
the plywood sheets with Px^’s at various intervals. The 
vertical members were placed four feet on centers in the 
shallow water section and from four feet to eighteen inches 
in the transition section. Additional rigidity was provided 
in the transition section by truss arrangements, as shown 
in Figure 4. Two 2x4 *s were attached to the top of the 
vertical studs and extended the length of the channel. 

A double floor constructed of two plywood sheets with 
two half-inch separators was used. The sides of the channel 
were bolted together, through the spacing of the double 
bottom, with three-eights inch threaded rod. 

All wooden sections exposed to water were water-proofed 
and painted with three coats of Morewear Vitri-Glaz 1320A. 

Dow Corning Sealant 780 was applied to all seams and joints. 

A paddle type wave generator was used to generate waves 
in the channel. An aluminum plate was hinged at the tank 
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Fig. 3 
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Fig.4 

TRANSITION SECTION TRUSS ARRG'T. 
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floor and attached by a pin connection at the top to a driving 
rod. Wave motion behind the paddle was damped by use of 
baffle plates. A two horse power variable speed drive, with 
an output speed range of twenty to one hundred eighty revolu- 
tions per minute, was mounted on top of the wave channel and 
fitted with a six and three quarter inch radius face plate. 
Ways were attached to the plate to adjust the driving rod 
eccentricity between zero and eight inches. 

The transition section contained a ramp with a slope of 
1:20 which shoaled the generated wave with a change in depth 
from four to two feet. The overall length of the transition 
section was 42 feet. 

The shallow water section had a depth of two feet and 
an overall length of 42 feet. The test module was located 
twenty four feet from the transition, allowing the waves to 
reach a fully developed state following the transition before 
reaching the test section. 

A variable slope dissipater beach was located at the end 
of the test section in order to dissipate the wave motion. 

The beach was constructed of metal shavings held between two 
pieces of perforated stainless steel sheet metal, which were 
separated by two inch wooden spacers. A solid sheet of one- 
eighth inch aluminum was attached to the bottom of the beach 
by three and five eighths inch separators. One by one inch 
aluminum angles were attached to the exposed surface of the 
beach, parallel to the wave fronts. Maximum slope of the 
beach was 1:7. 
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B . TEST MODULE 



The test module consisted of the circular cylinder with 
its own floor and walls, constructed as an integral system 
to provide easy removal of the model from the tank. The wave 
channel sides and floor were recessed to provide a smooth 
transition from channel to test module. The channel walls 
were made of plexiglass in the test module region to allow 
maximum visibility of the model. The test module is shown 
in Figure 5. 

A four inch diameter plexiglass cylinder was supported 
between the walls of the test module by use of cantilever 
beams and adjusted to approximately one-sixteenth inch above 
the floor. A thin flexible plastic barrier was installed in 
this gap to prevent any water motion under the cylinder due 
to wave motion. The barrier was held in place by ’O' ring 
material pressed into 0.007 inch slots in the cylinder and 
test module floor. 

The cantilever beams were used to support the cylinder, 
as shown in Figures 6 and 7. Bulkheads were fixed in the 
cylinder at approximately two inches from each end, and the 
fixed ends of the cantilever beams bolted on these bulkheads. 
The free ends of the two beams protruded slightly beyond the 
end of the cylinder and into the plexiglass test module walls 
where they were supported by small self-aligning ball bearings 
pressed into the plexiglass. Both beams were fitted with 
strain gauges and waterproofed using BLH Barrier ' C', as 
shown in Figure 8. 
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Pig. 5 TEST MODULE 
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Fig. 6 HORIZONTAL FORCE CANTILEVER 
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Pig. 7 VERTICAL FORCE CANTILEVER 
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Fig. 8 FORCE CANTILEVER 
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One beam, positioned with its largest cross sectional 
dimension vertical, was used to measure horizontal force and 
the other, with its largest cross section positioned hori- 
zontally, was used to measure vertical force. Calibration 
tests indicated that cross coupling was negligible. 

Dimensions of the beams were chosen to provide sufficient 
flexibility to allow measurement of forces using strain gauges 
and sufficient stiffness to provide a natural frequency that 
was large in comparison with the excitation frequency. In 
addition, the minimum width of the beams was constrained by 
the width necessary to mount the strain gauges. A reasonable 
trade-off was arrived at by making use of the maximum sensi- 
tivity of the amplifier/recorder and determining a maximum 
allowable value for the cantilever beam length. 

To prevent rotation of the cylinder, an offset arm was 
used to connect the horizontal force cantilever to the test 
module. A Farber Bearing, AMS 1K7, was used as a wheel at 
the test module end of the arm to reduce friction. 

C. WAVE HEIGHT PROBE 

• A capacitance type probe was chosen to measure wave 
heights in this study. The probe essentially consists of 
a single insulated wire which acted as one plate of a capac- 
itor, the water acting as the second plate. The capacitance 
varied as the water level rose and fell. A three-eighths 
inch diameter acrylic rod was attached to the bottom of the 
foil-like support member and the sensing wire was attached 
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at the bottom of the acrylic rod and insulated there. The 
tip of the acrylic rod was removable to provide easy changing 
of the sensing wire. The wire was connected at the top to 
a cable connector mounted on a nylon block. Number 30 A.W.G. 
wire with poly thermaleze insulation was used. 

The schematic of the electronic circuit utilized with 
this probe is shown in Figure 9. In this circuit, the 
square wave output of a multi- vibrator is used to drive a 
capacitance bridge, one leg of which is varied by the wave 
height probe. The resulting AC error signal is converted to 
a D.C. output signal by the full wave rectifier. The output 
signal was found to be insensitive to small multivibrator 
frequency changes and linear through the range of the wave 
height probe. A model LD 5211-13 Brush D.C. amplifier was 
used with the probe circuit. 

D. TEST PROCEDURE 

Calibration of the probe was conducted with the tank 
filled to the desired level and the test module and wave 
height probe in position. With the probe immersed to a depth 
of nine inches, the amplifier gain was set at desired sensi- 
tivity and the recorder pen bias adjusted to position the 
recorder to mid-scale. The probe immersion was varied over 
a seven inch range in half inch increments, using a traverse 
mechanism. This initial calibration was checked at several 
points before and after each set of runs. 
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Fig. 9 WAVE HEIGHT PROBE CIRCUIT 
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The force cantilevers were calibrated by loading the 
cylinder with a series of weights, using a pulley arrange- 
ment shown in Figure 10. The load was transmitted to the 
cylinder through a series of tapped holes located around the 
circumference at positions horizontally fore and aft, and 

vertical. The Brush recorder output was calibrated by loading 

\ 

the cylinder in half-pound increments. This calibration 
was also checked prior to and after each series of runs. 

After calibration was completed, the wave generator speed 
was set for the desired wave length and a data run commenced. 
During each run, wave generator eccentricity was varied from 
one inch to eight inches, or to a point where the waves 
broke in the channel. 
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FIGURE IQ. CANTILEVER CALIBRATION ARRANGEMENT 
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IV. PRESENTATION OF RESULTS AND CONCLUSIONS 



As discussed in the theoretical analysis, dimensionless 
wave force coefficients can be represented as functions of 
the dimensionless water depth (h/a), the dimensionless wave 

O 

height (H/2a), and dimensionless period (gT /h) parameters. 

A test program was developed which, for a given run, held 
water depth and period constant while varying wave height, 
in order to determine the effect of the dimensionless param- 
eters on the horizontal and vertical force coefficients. 

Four series of experimental runs were conducted at relative 
water depths of 9*0, 7.0, 5.5, 4.0. 

The dimensionless wave period parameter was varied for 
each series over a range of 100 to the maximum value attain- 
able with the experimental apparatus. The upper limit in 
wave height was usually provided by breaking. 

A series of representative traces are shown in Figures 
11, 12, 13 and 14 due to the difficulty in presenting all 
the traces. These traces represent experimental data for 
three dimensionless periods from each of the four dimension- 
less water depths. The uppermost trace in each run repre- 
sents wave height, the middle trace represents vertical 
force and the bottom trace represents horizontal force. The 
experimental values of the force coefficients for each 
experimental case are given in Appendix B. 

Using the dimensionless period, water depth and waveheight 
as inputs, the computer program was used to generate 

41 
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corresponding theoretical values of the force coefficients 
for each experimental case. This provided a comparison 
between experimental and theoretical values for each run. 
Plots of experimental and theoretical values for maximum 
horizontal and vertical force coefficients as a function of 
wave height for each period and water depth studied are 
included in Appendix C. 

Experimental results and theoretical computations 
indicate that, at the lower values of dimensionless period 
and wave height studied, the wave may be characterized as a 
linear wave. The horizontal force was found to vary sinu- 
soidally, and increased linearly with waveheight. The 
vertical force varied approximately as a negative cosine 
function in this region, having its maximum value in the wave 
trough. The magnitude of the vertical force coefficient 
was found to vary approximately as the waveheight squared. 
These characteristics are shown most clearly in runs Number 
44 through 47, 71 through 74 and 93 through 95 in Appendix C. 

As the relative period and wave height increased, the 
wave displayed the sharper peaks and longer, shallower trough 
that are characteristic of non-linear gravity waves. The 
horizontal force, proportional to the fluid particle acceler- 
ation, tended to reach maximum values closer to the wave 
peaks, but still remained symmetrical, with respect to the 
wave crest, as indicated by the theoretical traces. The 
theoretical analysis predicted a phase lag between the wave 
crest and the maximum force as small as 24 at the higher 



wave heights. 
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The vertical force component began to display positive 
values occurring at twice the frequency of the wave peaks, 
as shown in the experimental traces, Figures 11, 12, 13, 14. 

This is explained by the increasing importance of the lift 
component of vertical force, which, theoretically, is similar 
to a cosine squared function. At extremely long periods, the 
lift force is clearly dominant, and the maximum force no 
longer occurs in the middle of the wave trough. The maximum 
vertical force occurs when the cylinder is under the wave 
crests where the velocity is greatest. 

At large periods and wave heights, a point is reached 
where flow separation effects are no longer negligible. This 
is seen most clearly in the vertical force component, whose 
value becomes suddenly lower than the theory predicts, though 
still varying as the waveheight squared. This sudden drop 
in lift is apparently caused by a corresponding decrease in 
the lift coefficient resulting from flow separation. This 
phenomenon is most clearly evident in the graphical presentation 
of runs 51 through 56, and 71 through 73 in Appendix C. 

1). Horizontal Forces 

As previously mentioned, the horizontal forces 

tended to be nearly sinusoidal in nature with a maximum value 

which is quite linear with wave height. It is therefore 

possible to characterize the horizontal force coefficient 

2 

by it’s maximum value (F /pga £). The variation in the 

x max 

horizontal force coefficient with dimensionless wave height 
is shown graphically in Appendix C. 



The experimental values of horizontal force are seen 
to vary linearly with wave height for relative period param- 
eter values from 100 to 120. For relative period values 
greater than 120, the force coefficient varies in an increas- 
ingly non-linear fashion. This is consistent with the theory, 
as shown in Figure 15. However, for the range of period 
parameters studied, the deviation from linearity is relatively 
small. It is therefore possible to display the experimental 
data in terms of the slope of the force coefficient versus 
wave height plots, as shown in Figures 16 and 17. 

The theory predicts the maximum values of horizontal 
force coefficients prior to the wave heights where separation 
effects become apparent. The discrepancy between theoretical 
and experimental values is due to variation of the added mass 
coefficient in a separated flow as well as a net drag force 
which is no longer zero. 

2). Vertical Forces 

As indicated previously, the vertical force displays 
two regimes; one inertia dominated and one lift dominated. 

The maximum value of the vertical force coefficient 

p 

(F /pga 5.) is seen to vary as the wave height squared 
y max 

throughout the range of data taken. This is shown graphically 
in Appendix C. For higher values of relative wave height, 
the magnitude of the vertical force coefficients is less than 
predicted by the theory, due to separation effects, as 
explained above. 



48 



XO^OJO 




Fig. 15 HORIZ. WAVE FORCE 



49 




O O CD M" OJ 

• • • • 

CM — 



( D 2/H) 

(XVM) X j 



50 



10.0 20.0 40.0 60.0 100.0 200.0 400.0 

gT 2 /h 

Fig. 16 HORIZONTAL WAVE FORCE 




sz 

<v\ 

I — 
CD 



M/H IrPGj 



(xvw')\j 



51 



Fig. 17 HORIZONTAL WAVE FORCE 



B. CONCLUSIONS 



From the experimental and theoretical results, the 
following conclusions are considered warranted: 

1. The horizontal forces on the cylinder, resulting 
from the incident gravity waves, increased linearly with 
wave height for relative period values less than 120. For 
relative period values greater than 120, the force coeffi- 
cient was found to vary with wave height in an increasingly 
non-linear fashion. 

2. The vertical forces displayed two regimes; one inertia 
dominated at lower periods and wave heights, and one lift 
dominated at higher values of these parameters. 

3. For high period parameter values, separation effects 
can be expected for relative wave height values above 0.6, 
resulting in a sudden drop in experimental force coefficient 
values from those predicted by the theory. 

i 

4. The theory used in this study provides excellent 

2 

agreement with experimental values for relative period (gT /h) 
values from 100 to 200, and prior to the onset of flow 
separation. 
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APPENDIX A 



COEFFICIENTS FOR STOKES FIFTH ORDER WAVE THEORY 



S = sinh(d) 
C = cosh(d) 
tanh(d) 

( 8c 2| -8c 2 +9 ) 



(3840C 12 -4096C 10 +2592C 8 -1008C 6 +59 I t ! lC I| -l830C 2 +l47) 

512S 10 (6C 2 -1) 



1 

*isc 

(12C 8 +36C 6 -l62C 1< + l4lC 2 -27) 
192CS 9 



1 

S 

-C 2 (5C 2 +1) 

8S 5 

-(H84c 10 -1^^0C 8 -1992C 6 '+264lC lt -2^9C 2 +l8) 

1536S 11 

3 

8 ?*" 

(192C 8 -^24C 6 -312C 1 *+^80C 2 -17) 

768s 10 



33 



_ (13-4C 2 ) 
64S 7 



35 



^44 



( 512C 12 +4224C 10 -6800C 8 -12 ,8o8c 8 +l6 ,704c 11 - 3154C 2 +107) 

4096S 13 (6C 2 -1) 



(80C 6 -8l6C 1< -H338C 2 -197) 

1536S 10 (6C 2 -1) 



L _ -(2880C 10 -72 ,4 80C 8 +324 ,000C 6 -432, 0000^+163, 470C 2 -l6, 245) 
55 61 ,440S 11 (6C 2 -1) ( 8C -llC 2 +3 ) 



B 



22 



(2C +1) 
4S 3 



B 



24 



C(272C 8 -504C 6 -192C^+322C 2 +21) 



384S 



9 



B 



33 



3(8C°+1) 
6 3S 6 



B 



35 



B 



44 



( 88 , 128C l2j -208 , 224C 12 +70 , 8l8C 1Q +54 , 000C 8 -21 , 8l6C 8 +6264C^-54C 2 -8l ) 

12,228S 12 (6C 2 -1) 

C(768C 10 -448C 8 -48C 6 +48C 2J +106C 2 -21) 



B 



55 



384S 9 (6C 2 -1) 

(192,000C 18 -262,720C 1 ^+83,680C 12 +20,160C 1Q -7280C 8 ) + 
12,288S 19 (6C 2 -1) (8c 4 -iic 2 +3) 

+ (7l60C 6 -l800C 4 -1050C 2 +225) 

12,288S 10 (6C 2 -1) (8C^-llC 2 +3) 



A 1 = XA 1X + X \ 3 + X 5 A 15 
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A2 X ^22 ^ A 2 4 

A 3 = X 3 A 33 + x5a 35 
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DATA FOR RUN NUMBER 41. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 46. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 47. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 52. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 53. TOTAL NUMBER OF RUNS IS 23 
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DATA FQR RUN NUMBER 54. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 56. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 71. TOTAL NUMBER OF RUNS IS 23 
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OATA FOR RUN NUMBER720. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 73. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 74. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 75. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 92. TOTAL NUMBER OF RUNS IS 23 
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DATA FOR RUN NUMBER 93. TOTAL NUMBER OF RUNS IS 23 



oo 

Oo 

• o 

vO • 

H 

uoco 



x< 





(\! X 


UL 




• 


II 


UJ 




00 


* — 


0— 


»H >3* m 0 <M in co h O' >4- m c\j m 


UJ 


1-0 




OOOoOHHHMMfOvtinOvO 


X 


Om 


w 




O 


II 1- 


• 


000000000000600 


X 




UL 




►— 1 


oq: 






0 


0 






• 


MOO 






CO 


a:o 


U- 




pH 


UJl -4 


UJ 


OO'COCOinrHsOHfT^vOO'O'inCNJCO 




CLO 


0 


0 0 ^ h c\j co co -4- >t ^ <f* in vO 


CO 


< 


0 




♦—4 


_J<X 


— 


000060000000000 




< 


•X 




X 


20 


U-W 




1- 


01- 






CL 


»— 1 






tu 


cox 






0 


Zh 


< 


04 phh cor-co ino cor- 0 mo 0 




UJ CL 


CM 


f-4 c\i co co <t in I s - h- O' O' 0 c\j cm * 4 * in 


a: 


Xuj 


s 




UJ 


w O 


X 


OOoOOOOOOOhhhhh 


»- 


0 






< 








3: 










0 








X 


UJ — 


cocoHCOvoOvia^sOcocovf-a^mo 






o> 


0 0 0 H -h cm co co vt 0 r- 0 hco vf 






DL'w 




• 




O 


000000000 OO Ohh,-( 


00 


UJ 


LI- 




0 


a: 






X 


< 






0 








0 


CO 


LLI 


co in vO r^- co co 0 in vO 0 m 0 


UJ 


C£ 


O*-* 


ph »h cm co in so h- co 0 0 0 c\j >fr in o 


CO 


UJ 


O^X 




0 


1- 


Ow 


OOoOOOOoOHHiJHdH 


CO 


LU 


LL 






51 






• 


< 






CNJ 


or 








< 


h- 






CL 


X 


OCO^COhCOOoOOOOOOO 


CO 




LUO 


in co cm in O' <m O' o o O' co co o o o 




CO 


>•- 






CO 


< UJ 


oOfHHHc\ic\icococo*4-NMninO 


O 


UJ 


31 




O 


-J 






M 


X 








0 






UJ 


1—4 


1- 




CL 


CO 


zq: 


HCMcovtinvor^-coO^Or- (CMcOvtin 




X 


LULU 


O O O O O O O O O H r-4 tH H H H 


UJ 


UJ 


2 : CO 


cocococorococororocorocococofo 


> 


2 : 


OX 


O'O'O'O'O'O'O'O'O'O'O'O'O'O'O' 


< 


♦—4 


LUO 




3 


0 


CO ^ 





76 



DATA FOR RUN NUMBER 94. TOTAL NUMBER OF RUNS IS 23 
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APPENDIX C. WAVE HEIGHT-FORCE PLOTS 
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APPENDIX D 



COMPUTER PROGRAM 
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